Micro-evaporators for kinetic exploration of phase diagrams 
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We use pervaporation-based microfluidic devices to concentrate species in aqueous solutions with 
spatial and temporal control of the process. Using experiments and modelling, we quantitatively 
describe the advection-diffusion behavior of the concentration field of various solutions (electrolytes, 
colloids, etc) and demonstrate the potential of these devices as universal tools for the kinetic explo- 
ration of the phases and textures that form upon concentration. 
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Determination of the phase diagram of multicompo- 
nent systems is of importance in many realms: industrial 
formulation, protein cristallization, bottom up material 
assembly from spontaneous ordering of surfactant, poly- 
meric or colloidal systems [H 0- Q Depending on the 
application, one may want to access only the equilibrium 
phase diagram or gain additional information as to the 
metastable phases that can appear for kinetic reasons. 
Methods to reach these goals often imply tedious and sys- 
tematic measurements, requiring for screening purposes 
the use of robotic platforms. Two generic strategies con- 
sist in varying (in space or time) the temperature of sam- 
ples of given concentrations on the one hand, and on the 
other hand isothermal concentration by cither removal 
of the solvent (osmosis, drying), external action on the 
solutes (sedimentation or dielectrophorcsis for colloids), 
or studies of spontaneous intcrdiffusion in contact exper- 
iments. 

In this Letter we introduce microfluidic tools for con- 
trolled isothermal concentration of a wide range of sys- 
tems, covering solutions of ions, polymers, proteins, sur- 
factants and colloidal suspensions. Our work is inspired 
by recent observations 0, 0| that in standard microsys- 
tems built of PolyDiMethylSiloxane (PDMS), sponta- 
neous water permeation through the PDMS matrix in- 
duces flows that can be used to concentrate colloids. Tak- 
ing a step further, we have engineered specialized micro- 
geometries that allow us to control spatially and tem- 
porally the evaporation process as well as the resulting 
concentration of solutes. Their parallel implementation 
in microfluidic format could open the way to fast screen- 
ing methods. 

After a brief description of the micro-devices, we 
demonstrate first our control of the concentration process 
on dilute aqueous solutions of fluorescein and nanopar- 
ticles in a simple geometry. We then discuss how mi- 
crofabrication permits to widen the range of possibili- 
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FIG. 1: Sketch of the finger geometry: top and side views 
showing the gas and liquid layers and the thin PDMS mem- 
brane in between (typical dimensions: e = 10 /im, h = 20 (im, 
w = 200 fim, Ln = 10 mm). 



ties and applications of such devices. As a study case, 
we report controlled nucleation and growth of crystals 
of potassium chloride (KC1), and show how such experi- 
ments provide quantitative information on various ther- 
modynamic quantities (solubility, crystal density) as well 
as kinetic features (sensitive to the rate of concentration) . 

The devices - The devices used in this paper are two- 
layer PDMS on glass microsystems (Fig. ^) (fabrica- 
tion procedure detailed in jjj). The microchannels of 
the bottom layer are filled with the solution of interest, 
while air (at controlled humidity) is circulated through 
the microchannels of the top layer so as to remove the 
water that pervaporates through the thin membrane of 
PDMS that separates the two networks where they over- 
lap (thickness e in the 10 — 30 \xm range). Many com- 
binations of geometries for the bottom and top networks 
can be envisaged. We focus here on the simple "finger" 
geometry of Fig. 1, a dead end channel of rectangular 
cross section (height h, width w, length L) connected to 
a larger (millimetric) feeding reservoir containing the so- 
lution to be concentrated. A terminal section of length 
Lq < L (typically mms to cms) is covered by the water 
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FIG. 2: Velocity of fluorescent tracers at a fixed location in 
a finger in response to the switching on and off of air circu- 
lation in the water removal network. Instantaneous velocity 
(gray lines) is obtained from individual trajectories gathered 
by Particle Tracking Velocimetry on 1.1 /im diameter trac- 
ers. The symbols are for the mean velocity averaged on ~ 10 
trajectories. Lo = 12.5 mm, h — 22 ^m, e ~ 20 /im. A la- 
tency time of « 0.2 s is observed before measurements during 
which the PDMS membrane undergoes vibrations due to the 
pressure switch. 



removal network. 

The operation principle is simple: water in the bottom 
channel pervaporates through the thin membrane, which 
induces a compensating flow from the reservoir and con- 
centration of the convected solutes at the tip of the finger 
0,0. This is similar to concentration at the boundary of 
a drying droplet 0, , without the motion and shear of 
the concentration zone due to the recess of the liquid-air 
interface, and without spurious convective flows, thanks 
to the confinement. The small dimensions lead to ther- 
mal regulation much faster that typicall kinetics involved 
in our studies (ms vs hours), which permits isothermal 
studies. 

Control of flow and particle concentration - To quantify 
the induced flow, we adapt the analysis of 0, , antici- 
pating that evaporation occurs here mostly through the 
thin membrane, at a volumic flow rate of water v e (it 
has dimension of a velocity). Mass conservation then 
sets the height-averaged velocity in the micro channel: 
v(x) = —v e (x/h). Its amplitude rises with the distance 
x from the dead-end (x = 0) up to vq = v e (Lo/h) at the 
end of the evaporation zone (x = Lo)- v(x) = —vq in the 
evaporation- free section of the finger Lq < x < L. 

Compared to previous work 0, , with our dedicated 
systems we induce larger values of v e (in the 50 nm/s 
range) due to the thin membranes (permeation yields a 
limit scaling as 1/e) and to the dry air flown through 
the top network (velocities of order cm/s) that reduces 
the diffusive boundary layer and maintains constant the 
driving force fpr evaporation. More importantly, we gain 
a spatial and temporal control on v e by designing the 
geometry (evaporation is negligible but for chosen loca- 
tions) and by tuning in time the air flow and thus v e . 
Quantitative temporal control of the flow field is clear 
from the motion of tracers at a given location as the air 



flow is successively turned on and off (Fig. 2). When on, 
tracer velocities of order vq = 13 /im/s are observed, cor- 
responding to r e = h/v e ~ 10 3 8 and v e ~ 22 nm/s. The 
velocity drops below 1 pm/s when the air flow is turned 
off, after a response time of a few seconds, compatible 
with that of the water flux through the thin PDMS layer, 
e 2 /DpDMS ~ 0.5 s for e ~ 20 /im and a diffusion coeffi- 
cient for water in PDMS £>pdms ~ 10~ 9 m 2 /s d Gil- 
Control of the flow field translates into that of the in- 
duced concentration process. For the simplest case of 
a dilute species of diffusion coefficient D in a finger, in 
a one-dimensional description the conservation equation 
dtc + d x J = relates the concentration c(x,i) and flux 
J(x, t) = cv — Dd x c. We focus now on steady evaporation 
and thus constant v(x) in time, with a reservoir at fixed 
concentration cq- The physics at work is simple: the flow 
converts the solute towards the dead end where it accu- 
mulates. The current of solute injected into the finger is 
steady Jq = cqvq. Backwards thermal diffusion against 
the flow controls the width of the accumulation zone [f| , 
which is p — (Dh/v e y/ 2 = {Dt^ 1 / 2 < L for strong 
flows or long fingers vqL /D 3> 1 (which we assume from 
now on |l lj). At distances larger than p, diffusion is negli- 
gible and the suspension is simply concentrated by water 
removal at constant particle flux c(x)v(x) = cqVq. Al- 
together, after a transient of duration p 2 /D w r e , the 
profile is well approximated by a Gaussian (because of 
the linearity of v(x)) increasing linearly in time, "fed" by 
a steady hyperbolic ramp delivering a current Jo = cqVo- 
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with R(x) ~ Lq/x for p <C x <C Lq. Experiments on 
solutions of fluorescein and nanocolloids quantitatively 
support this analysis (Fig.|2J). 

We thus have multiple ways to reach quantitative 
control of the kinetics of the process. Indeed, the 
concentration at the tip increases as 4f(x = 0) = 
T i /2/n(v e /h) 3 / 2 CoLoD~ 1 / 2 , which shows that control can 
be exerted through geometrical features Lq, h and e (that 
affects v e ), and through operational parameters Co and 
v e . The latter can be modulated during an experiment, 
and this allows us to pinch or relax concentration pro- 
files at a rate which we have checked to be the diffusion- 
limited response time ~ p 2 /D = h/v e = r e (independent 
of the species), much larger than the response time of the 
flow (typically minutes instead of seconds). 

Benefits from micro fabrication - Beyond the simple 
PDMS finger discussed above, microfabrication offers 
many possibilities that we are currently exploring : 
- We can parallelize tests so as to perform rapid screen- 
ing with minute amount of material. Indeed numerous 
microchannels of various geometries and types can be 
fabricated on a single chip, and can be fed from a single 
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FIG. 3: Evaporation-induced concentration. Top: log-linear 
plot of fluorescence intensity against position at different 
times for a finger filled with a aqueous solution of fluorescein. 
The fit corresponds to the prediction Q. The concentra- 
tion in the accumulation zone is well described by a Gaussian 
hump of width p ~ 500 /im, that increases linearly in time 
(insert), yielding v e — 15 ± 2nm/s and r e = 950 ± 100 s. 
(co = 5 10 M, Lo = 4.7 mm, h = 15 fim, e ~ 25 /im; data 
corrected from photobleaching with a characteristic time of 
~ 2400 s and total exposure time of « 60 s). Bottom: images 
from four fingers in similar conditions showing that the size 
of the accumulation zone varies with the diffusion coefficient 
of the markers. 



reservoir or multiple ones. We touch upon the corre- 
sponding potentialities in the experiment reported below 
where a 15-finger chip is used (Fig. HJ. 
- Geometries are not limited to fingers. Let us illustrate 
how this can help on an important example. Equation (1) 
indicates a drawback of the finger geometry for the anal- 
ysis of ill-characterized mixtures (a drawback actually 
shared by many other methods : drying, field induced 
concentration, ...): different solutes are concentrated at 
different rates {p depends on D) so their proportions in 
the accumulation zone differ (in an a priori unknown way) 
from those in the reservoir. This can be alleviated by at 
least two solutions relying on appropriate geometrical de- 
sign of the channels. The first is to fabricate very long 
serpentine fingers and focus on the hyperbolic concen- 
tration ramp (R(x) in equation (1)) where all species are 
concentrated alike. A second strategy consists in adding 
a chamber of area A (thickness h) at the end of the finger 
with pcrvaporation limited to the finger. The area of the 
accumulation zone is then ~ A+pw, and thus essentially 
independent of p (or D) as long as p <C A/w (with A a 
few mm 2 and w a few hundred microns, this encompasses 
anything from ions to micron-sized colloids). Although 



the concentration rate is lower than in a linear geometry 
due to the increased accumulation area, the various so- 
lutes concentrate at the same rate. 

- The PDMS devices presented here will work only with 
a limited set of solvents that do not swell this elastomer. 
However microfabrication permits the extension of the 
concepts presented here to other materials. The bottom 
layer of channels can for example be etched in an im- 
permeable solid matrix (e.g. glass for the "floor" and 
"side- walls" ) separated from the air-circulating network 
by a thin membrane of a chosen porous materials. Similar 
sandwich micro-constructions have been used in related 
contexts HlHIll. 

Controlled crystal nucleation and growth - Let us now 
illustrate how such systems can be used to study physical 
phenomena that occur upon concentration, by focusing 
on a well characterized system, KC1 aqueous solutions. 
We use a microchip with multiple "fingers" of different 
lengths Lo originating from the same reservoir, and per- 
form a few experiments, at the same steady air flux, with 
different intial concentrations. Upon concentration, we 
observe in each finger the nucleation of crystals close 
to tip, at x c and t c (Fig. 0J, and then their subsequent 
growth with a front location at Xf(t). The time scales 
involved vary widely with initial concentration en and fin- 
ger length Lq (Fig. BJ. Nevertheless, a rescaling of data 
for both nucleation and growth is possible if we account 
for the concentration rate proportional to cqLq as we now 
show. 

Figure 21 shows that nucleation occurs at a typical 
time that scales t c = K^Lq) -1 , with K a constant. 
From (1) this corresponds to a nucleation concentra- 
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tion c c = Ky/zJ^D- 1 ' 2 ^ 
KC1 (c c = 4.0 ± 0.5 M, D ~ 2.0 ± 0.2 10~ 5 cm 2 /s at 
this temperature) are consistent with the experimental 
value of K ~ 5.5 ± 1 sMm and a reasonable value for 
r e = h/v e = 800 ± 50 s (which yields p w 1200 /im). The 
rescaling of crystal growth data in Figure 0] (middle), 
suggests that the process is limited by the solute feed at 
Jo = cqVq = C()Lo/r e . Indeed in such a picture, if we as- 
sume that the growing crystal at solute concentration cx 
fills the channel, then the initial growth rate is dictated 
by mass conservation —J- = Jq/cx = (cqLq)/ (r e cx). 

The observed value (coio)" 1 ^ 1 = 0.6 ± 0.1 10~ 4 sMm 
yields a value cx = 20±2 M consistent with the literature 
value (26 M). The observed decrease of the growth rate 
in time for each channel likely results from that of the 
evaporation length as the crystal invades the channels. 

The above analysis demonstrates a possible use of our 
devices. Feeding a few fingers with a system of reference 
(e.g. KC1 or fluorescent probes) provides a local calibra- 
tion of the value of r e (to account for variations from 
device to device). Then measurements of t c and front 
motion in other fingers fed with the solution to analyze 
yield estimates of c c D x / 2 and of the solute concentration 
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FIG. 4: Top: view of the chip with 15 micro channels with 
different evaporation lengths Lo (the bottom of the frame 
is the end of the evaporation zone). In each channel, the 
evaporation-induced concentration leads to crystalization of 
KC1 solutions at x c after a time t c . Crystal growth is then 
monitored Xf (t). Middle: Left: crystal fronts Xf(t) for various 
evaporation lengths Lo and initial concentrations Co = 37 mM 
(bottom), c = 123 mM (middle), c = 372 mM (top). Right: 
plots of nucleation time t c against 1/coLo, and front growth 
Xf against coLo(t — t c ). Bottom: Crystals growing at various 
Co: "droplets" in front of the crystal are only visible at low 
concentrations. 

cx in the phase that nucleates. 

Further, with our micro-devices we can investigate the 
influence of the kinetics of the concentration process on 
the morphologies and phases produced. We indeed ob- 
serve on KC1 solutions that both the nucleation and 
growth processes qualitatively change when the initial 
concentration Co (and thus the concentration rate) is var- 
ied |llj . In the studies reported above, we reproducibly 
observe what looks like liquid droplets in the dense elec- 
trolyte solution in front of the growing crystal, but only 
at the lowest initial concentration examined (Fig. ^bot- 
tom) , a phenomenon that we are currently investigating. 
Remarkably, the various transient organizations (nucle- 
ation scenario, presence or absence of droplets) do not 
affect the rescalings of Fig. 4, possibly due to the nar- 
row metastable region of KC1 and the robustness of mass 
conservation arguments. 

Perspectives - Although we have used KC1 to demon- 
strate some of the potentialities of such microevapora- 
tors, their main interest likely lies in the rich world of 



soft matter, where systems adopt many diverse and less 
known forms of organization |l5l Hf3 | . Beyond ionic so- 
lutions and colloids, we are currently performing experi- 
ments on surfactant solutions and on proteins Our 
early experiments on solutions of Sodium Dodecyl Sul- 
fate (SDS) indeed allow us to observe up to four different 
mesophases in a given microchannel, and more impor- 
tantly complex and intriguing dependences of that se- 
quence and of the texture of each phase on the kinetics 
and history of the concentration process. At the same 
time, data for nucleation and front growth obey rescal- 
ings akin to those of KC1, showing our kinetic control. 

For all the reasons presented above, and especially the 
delicate control on the kinetics and their wide applicabil- 
ity to mixtures of solutes anywhere between simple ions 
to micron-sized colloids, we believe that our microfluidic 
evaporators are powerful and rather universal means to 
explore the wonders of soft matter. 
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